On the microcanonical solution of a system of fully coupled particles 
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We study the Hamiltonian Mean Field (HMF) model, a 
system of fully coupled particles, in the microcanonical en- 
semble. We use the previously obtained free energy in the 
canonical ensemble to derive entropy as a function of en- 
ergy, using Legendre transform techniques. The temperature- 
energy relation is found to coincide with the one obtained 
in the canonical ensemble and includes a metastable branch 
which represents spatially homogeneous states below the crit- 
ical energy. "Water bag" states, with removed tails momen- 
tum distribution, lying on this branch, are shown to relax to 
equilibrium on a time which diverges linearly with N in an 
energy region just below the phase transition. 

PACS numbers:05.40+j,05.70Fh 



I. INTRODUCTION 

In a pioneering work Fermi, Pasta and Ulam jlj showed 
that in a Hamiltonian system with short-range interac- 
tions the relaxation time to thermodynamical equilib- 
rium can be extremely long if the energy is small enough. 
This striking feature triggered an intense research ac- 
tivity and, nowadays, the relation between macroscopic 
thermodynamical properties and microscopic dynamics 
is actively investigated in a large variety of physical sys- 
tems g. The common belief that microcanonical time 
averages coincide with equilibrium ensemble averages af- 
ter a possibly quite long transient is confirmed by sev- 
eral numerical experiments (see Refs. || for reviews). 
Therefore, although the rigorous proof of ergodicity is re- 
stricted to a few rather abstract dynamical systems (e.g. 
the Sinai billiard), it is nevertheless believed, and con- 
firmed numerically 0j5|, that time averages of thermo- 
dynamic functions (temperature, internal energy, specific 
heat) converge to their equilibrium values, at least in the 
regime of strong chaos. 

The case of long-range forces is quite less explored from 
this point of view. The study of relaxation to thermal 
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equilibrium is here made more complex by the conse- 
quences of the ineguivalence of statistical ensembles. In- 
deed, it is well known that the impossibility to isolate a 
small sub-system from a large thermal bath, due to the 
presence of long-range forces, prevents the derivation of 
canonical from microcanonical ensemble. 

Exactly solvable toy models (6) with weak long-range 
forces explicitly display ensemble inequivalence and the 
presence of negative specific heat in the microcanonical 
ensemble. Similar phenomena were also obtained for 
gravitational models 0. These toy models show a first 
order phase-transition in the canonical ensemble, which 
corresponds to the presence of a concave region of the 
entropy-energy curve in the microcanonical ensemble. 
This feature is present also in microcanonical molecu- 
lar dynamics simulations of short-range models without 
hard-core JsJ , where one observes, in the phase-transition 
region, a high-density clustered phase coexisting with a 
low-density gaseous phase. More recently, microcanoni- 
cal Monte-Carlo simulations have shown similar proper- 
ties for the diluted three-state Potts model Q. 

All the previous remarks concern equilibrium proper- 
ties, but it is known that long-range interacting systems 
show also extremely slow relaxation to equilibrium. A re- 
markable example is the one-dimensional self-gravitating 
potential, where the study of long transients and quasi- 
cquilibria has been undertaken by many groups p0[ . 
More recently, a scaling law for the relaxation to equilib- 
rium has been proposed for gravitational systems based 
on diffusion in phase-space . 

The Hamiltonian Mean Field (HMF) model @ has 
been introduced with the aim of studying clustering phe- 
nomena in A/-body systems in one dimension. The origi- 
nal motivation was to consider a truncated Fourier series 
development of the one-dimensional gravitational and 
charged sheet models. Its canonical ensemble exact so- 
lution predicts a second-order phase transition from a 
clustered phase to a gaseous one. HMF is an infinite 
range interacting system and one would naturally ex- 
pect unusual thermodynamical behaviour. However, mi- 
crocanonical simulations of HMF show that "carefully" 
prepared initial states lead to a good agreement with 
canonical analytical results for the energy-temperature 
relation On the other hand, "water bag" initial 

states, ie. states with momentum distribution on a fi- 
nite support, are metastable and show an extremely slow 
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relaxation to equilibrium near the critical energy p3| . 

In this paper, we derive analytically the entropy-energy 
relation of the HMF model in the microcanonical ensem- 
ble. We further study the properties of metastable states 
corresponding to the low energy extension of the high en- 
ergy gaseous phase. The relaxation time of these states to 
canonical equilibrium is also investigated, showing that 
it increases linearly with the number of particles N near 
the critical energy. 



II. THE MODEL 

The Hamiltonian Mean Field model describes the dy- 
namics of N classical and identical particles confined to 
move on the unit circle [^2| . Its Hamiltonian writes 

N 2 N 

F = Ey + ^E( 1 - C0S ^-^))=^ + y > C 1 ) 

where (9i,pi) are canonically conjugated coordinate and 
momentum variables with 9i 6 [0, 1-k\. K is the kinetic 
energy and the potential energy V corresponds to the first 
harmonic of the one-dimensional self-gravitating poten- 
tial V = J2i \8i ~ @j\ ( U P to scaling constants). Total 
momentum 22 i pi is also conserved and always set to zero 
without loss of generality. 

In the canonical ensemble, model (0) undergoes a 
second-order phase transition at the inverse temperature 
value p c = T~ l = 2 fl2fl , corresponding to a critical spe- 
cific energy U c — E c /N = 3/4. The free energy 
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where 



Z(j3,N) = jY[d Pi deiexp(-PH) 
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has been obtained using saddle-point techniques flj2| , giv- 
ing 
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whose unique solution at small j3 is x = 0, while at large 
(3, ie. low temperature, a pair of symmetric non vanishing 
solutions is present. The bifurcation (of pitchfork type) 
occurs at f3 c = 2 and at this temperature the value of 
the free energy F has a discontinuity in the second order 
derivative, which is the signature of a second-order phase 
transition. This is also confirmed by the calculation of 
the order parameter 
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which vanishes continuously at /3 C) remaining zero in the 
whole high temperature phase. At low temperatures, M 
measures the degree of clustering of the particles. 

The equations of motion derived from model (Q) , writ- 
ten in terms of M, show that each particle has a pendu- 
lum like motion 
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which makes explicit the self-consistent nature of the dy- 
namics, since each particle moves in the M-field, which 
is itself determined by the position 6>; of all the parti- 
cles. This has led [[l2| to the interpretation of the phase 
transition as a dynamical process of particle evaporation 
from the cluster. Moreover, this approach follows the 
idea that the study of self-consistent TV-body dynamics 
is simplified by an "effective" reduction of the number of 
degrees of freedom 



III. ENTROPY - ENERGY RELATION. 

The free energy F{0) in the canonical ensemble is read- 
ily computed once one knows the dependence of entropy 
S on specific energy U = E/N in the microcanonical 
ensemble. One begins from the relation 



Z{P,N)= / uj(E,N) e - (3E dE, 



where 



u(E, N)=[Yl dpid9 t 5(E - H) 



(9) 
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where 



X{x,(3)=\n{2-Kl {x))- — 
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7o being the zero order modified Bessel function. The 
maximization of X is performed by solving the consis- 
tency equation 



x I\ (x) 

= 



(6) 



is the microcanonical phase-space density at energy E. 
The lower limit of the integral in (^) is given by the energy 
of the ground state, which vanishes in our model (in this 
case all of the particles are at the same position with zero 
momentum). After rescaling by N and exponentiating u 
Eq. (0) becomes 



Z(0,N)=N I dU exp(N(-(3U + ^\n(uj(E,N))) 



(11) 
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This integral is then solved by the saddle-point technique 
in the N — > oo limit, which implies taking the maximum 
of the argument of the exponential. Using the definition 
of entropy 



S(U) 
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and recalling the definition of free energy (g), one finally 
gets 



(3F(P) = max[-/3C/ + S(U)} 



(13) 



Entropy S and free energy F are thus related by a Leg- 
endre transform. Assuming now that S is concave (i.e. 
downward bended)]^], one can invert (|l3|) and obtain 



S(U)=mm[f3(U-F((3))}- 



(14) 



This assumption is not innocent and the invertion 
( ETjj ) can be safely used only for second order phase- 
transitions, as it is for our model. On the contrary, it 
is well known that for first-order phase transitions en- 
tropy has a concave region |^,^| , which is responsible for 
the negative specific heat in the microcanonical ensemble. 
A generalization of model (|l|) to particle motion on the 
two-dimensional torus |L5) does have a first-order phase 
transition, making the invertion impossible in this case. 

The calculation of entropy in the microcanonical en- 
semble is then reduced to the min-max procedure 



S(U) = minmax[/3C7 + - ln(^) 

/3>0 x 2 P 
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-X(x,{3)] , (15) 



with X given by @ . For (3 > (3 C = 2 the max is at x = 
and the min at (3 = 1/ (2U — 1), which gives the entropy 
function 



S(U) = - ln(2[7 - 1) + const. 



(16) 



This is plotted in fig. [l] with full line above U = U c = 3/4 
and with dashed line for 1/2 < U < U c (below U = 
1/2 the temperature T — 1/(3 becomes negative and the 
solution meaningless). This solution corresponds to the 
high temperature Homogeneous Phase (HP) with M = 
for U > U c and to the MetaStable (the meaning is 
clarified below) Homogeneous Phase MS hp below U c . 
Both phases are present also in the canonical solution 
and while in this latter context the MShp corresponds 
to a local, but not global, minimum of the free energy, 
in the microcanonical context the entropy of the MShp 
is not maximal. Indeed, the solution which maximizes 
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the entropy bifurcates continuously from this one at U c . 
Below this energy the equilibrium solution is related to 
the two symmetric solutions with x ^ of the consistency 
equation (j^). They give the new maximum value of X 
in ( |l5| ) and the numerical search of the minimum in [3, 
which is unique, allows to compute the stable branch of 
S(U) below U c . The latter is plotted in fig. |l| with the 
continuous line and labelled CP, because it corresponds 
to the Clustered Phase, which is characterized by a non 
vanishing value of M. The temperature T can also be 
derived from the standard relation dS/dU = T"" 1 and is 
reported in fig. The stable solution S(U) (full line) 
is continuous with its derivative. The discontinuity is in 
the second derivative, which gives a jump in the specific 
heat C v = (dT/dU)- 1 . 

To summarize, our results suggest that the microcanon- 
ical solution is fully consistent with the one in the canon- 
ical ensemble, and we can therefore conclude that the two 
ensembles lead to equivalent results for model nA). 



IV. METASTABILITY 

However, the dynamics reveals a quite peculiar behav- 
ior of the microcanonical metastable branch. We have 
performed numerical simulations in microcanonical fixed- 
energy conditions by integrating the equations of motion 
(||), starting from appropriately chosen initial states. We 
have used a fourth order symplectic algorithm fl6| with 
a time step 0.1, allowing relative energy conservation as 
good as 10 -5 . We have checked the energy-temperature 
curve, definining temperature as twice the time average 
of the kinetic energy: T = 2 < K >. The equilibrium 
energy-temperature solution is reached from generic ini- 
tial states over most part of the energy range. However, 
for 1/2 < U < U c relaxation to equilibrium becomes very 
slow. We have started the runs in this range with the par- 
ticles homogeneously distributed on the circle (M = 0) 
and with a "water bag" momentum distribution: the mo- 
menta are uniformly distributed in an interval around 
zero such that the variance of the distribution coincides 
initially with the temperature of the metastable branch 
T = 2U-1. 

The relaxation of the initial "water bag" to the equi- 
librium state proceeds through the slow development of 
Maxwellian tails in the velocity distribution. In the lit- 
erature one can already find indications that the relax- 
ation time to equilibrium grows linearly with N jl3| (the 
authors used as a criterion the convergence to the Boltz- 
mann entropy of the Maxwellian). We have decided to 
consider directly the time evolution of the temperature 
of the MShp states. The initial temperature lies below 
the equilibrium one (see dashed line in fig. [l]) and re- 
laxes to a fixed threshold temperature o~T eq after a time 
~tre.ia.x-, where a < 1 is the threshold amplitude and T eq is 
the microcanonical equilibrium temperature represented 
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by the full line in fig. [|. Considering different initial 
states, we have observed that the fluctuations of t re i ax 
are Poissonian and that a well defined average exists. 
This average relaxation time t av _ re iax is plotted in fig. ^ 
for a = 0.9 as a function of iV for different values of 
U. For U < 0.56, t av _ re i ax is almost independent of N 
whereas for larger values of U, t av _ re i ax increases approx- 
imately as N. Nothing substantially changes using dif- 
ferent thresholds by changing the value of a. We are 
therefore able to confirm, for a different observable and 
a better statistics (we have considered up to 3600 initial 
conditions for the points in fig. ||), the results of Ref. |b| 
on the linear divergence of the relaxation time with N, 
at least in a restricted energy interval just below U c . The 
existence of such divergence implies that performing the 
N — > oo limit before the t — > oo limit, the system will 
stay indefinitely out of equilibrium and the momentum 
distribution will never develop Maxwellian tails. This has 
attracted the attention of Tsallis fTjJ , who argues that his 
entropy definition should be appropriate to describe the 
thermodynamical properties of such metastable states in 
the thermodynamic limit (in particular it would possibly 
describe the power law tails of the momentum distribu- 
tion) . It should however be remarked that what we have 
tested here is not the divergence of the lifetime of such 
states, but that of the relaxation time. However, recent 
numerical simulations performed by Rapisarda pl| seem 
to indicate also the divergence of the lifetime with N, en- 
couraging Tsallis' interpretation. Furthermore, it is pos- 
sible to prepare states on the metastable branch which 
have initially a Maxwellian momentum distribution and 
those states are observed to relax to equilibrium pretty 
fast. It seems that the "water bag" condition, which 
corresponds to an initial truncation of the tails, is a nec- 
essary one for obtaining a divergence of the relaxation 
time. 

V. CONCLUSIONS AND PERSPECTIVES 

Let us summarize what we have found. The micro- 
canonical solution of the Hamiltonian Mean Field model 
is here obtained with Legendre transform techniques, 
which avoid state counting. The entropy is computed 
using a min-max procedure and the energy-temperature 
relation is derived analytically. The equilibrium states 
are shown to coincide with the canonical ones in the 
thermodynamic limit. This result is not obvious since, 
due to the long-range intcrparticle coupling, ensemble 
equivalence might have been violated. Besides the equi- 
librium solutions, other states are found, corresponding 
to the metastable homogeneous configurations lying be- 
low the critical energy. The study of their relaxation 
time to equilibrium exihibits a divergence proportional 
to the number of particles, meaning that in the thermo- 
dynamic limit the trajectory in the phase-space remains 



indefinitely trapped in a non-equilibrium region (some- 
thing similar to so-called "mechanical" initial states for 
the Fermi-Pasta- Ulam model in Ref. [|l9|). Further in- 
vestigations are necessary to understand the relevance of 
the metastable states for other more realistic long-range 
forces, beginning with the very interesting extension of 
the HMF proposed by Tsallis and co-workers p0| , pl[ , 
whose canonical solution has been recently obtained in 
Ref. Q . These results could be relevant for the physics 
of self-gravitating systems and for nuclear dynamics, 
where long-range forces play a distinctive role. 



Acknowledgements 

We thank V. Latora, A. Rapisarda and A. Torcini for 
fruitful discussions and ongoing collaboration on this sub- 
ject. M. A. acknowledges D. Noack for important logis- 
tics support. MA. and S.R. thank the Institute for Sci- 
entific Interchange, Torino (Italy) and the Laboratoire de 
Physique of the ENS-Lyon (France) for hospitality and 
financial support. The Max-Planck Institute for Physics 
of Complex Systems in Dresden (Germany) is also ac- 
knowledged for allowing us to use its powerful computer 
center. 



[1] E. Fermi, J. Pasta and S. Ulam, Los Alamos Report LA- 

1940 (1955), in D.C. Mattis, The Many-Body Problem, 

World Scientific, (1993). 
[2] Chaos, Focus issue on Chaos and Irreversibility, T. Tel, 

P. Gaspard and G. Nicolis, Eds. (1998). 
[3] S. RufTo, in Symmetry and perturbation theory, D. Bam- 

busi and G. Gaeta eds., GNFM-CNR (1998); P. Poggi 

and S. Ruffo, Physica D, 103, 251 (1997). 
[4] R. Livi, M. Pettini, S. Ruffo, A. Vulpiani, J. Stat. Phys., 

48, 539 (1987); D. Escande, H. Kantz, R. Livi and S. 

Ruffo, J. Stat. Phys., 76, 605 (1994). 
[5] C. Giardina and R. Livi, J. Stat. Phys., 91, 1027 (1998). 
[6] P. Hertel and W. Thirring, Ann. of Physics, 63, 520 

(1970). 

[7] D. Lynden-Bell and R. Wood, Mon. Not. R. Astr. Soc, 

138, 495 (1968). 
[8] A. Compagner, C. Bruin and A. Roelse, Phys. Rev. 

A, 39, 5989 (1989); HA. Posch, H. Narnhofer and W. 

Thirring, Phys. Rev. A, 42, 1880 (1990). 
[9] D.H.E. Gross, A. Ecker and X.Z. Zhang, Ann. Physik, 

5, 446 (1996) and contribution of D.H.E. Gross to this 

volume. 

[10] H.L. Wright, B.N. Miller and W.E. Stein, Astrophys. 
Space Sci. 84, 421 (1982); B.N. Miller and C.J. Reidl 
Jr., Astrophys. J. 348, 203 (1990); T. Tsuchiya, N. 



4 



[11] 

[12] 

[13] 
[14] 

[15] 

[16] 
[17] 

[18] 
[19] 

[20] 

[21] 

M 



Gouda and T. Konishi, Phys. Rev. E, 53, 2210 (1996); 
T. Tsuchiya, N. Gouda and T. Konishi, Astrophys. Space 
ScL, 257, 319 (1997). 

M.-L. Chabanol, F. Corson and Y. Pomeau, Europhys. 
Lett. 50, 148 (2000). 

S. Ruffo, in Transport and Plasma Physics, S. Benkadda, 
Y. Elskens and F. Doveil eds., World Scientific, Singa- 
pore, 114 (1994); M. Antoni and S. Ruffo, Phys. Rev. E 
52, 2361, (1995). 

V. Latora, A. Rapisarda and S. Ruffo, Phys. Rev. Lett., 
80, 692 (1998); Physica D, 131, 38 (1999). 
J.L. Tennyson, J.D. Meiss and P.J. Morrison, Physica D 
71, 1 (1994); D. del-Castillo-Negrete, Phys. Lett. A 241, 
99 (1998). 

M. Antoni and A. Torcini, Phys. Rev. E 57, R6233 

(1998) , A. Torcini and M. Antoni, Phys. Rev. E 59, 

(1999) . 

I. McLachlan and P. Atela, Nonlinear ity 5, 541 (1992) 



C. Tsallis, Braz. J. Phys. 29, 1 (1999) , |cond-mat /9903356 
and refs. therein. 

A. Rapisarda, contribution to this conference. 

J. de Luca, A.J. Lichtenberg and S. Ruffo, Phys. Rev. E, 

60, 3781 (1999) 

C. Anteneodo and C. Tsallis, Phys. Rev. Lett. 80, 5313 
(1998). 

F. Tamarit and C. Anteneodo, Phys. Rev. Lett. 84, 208 
(2000). 



A. Campa, A . Giansanti and D. Moroni, preprint cond 
mat/0002168 and contribution of A. Giansanti to this 



* * U=0.54 






o o XJ=0.56 






♦ ♦ U=0.5S 


~ 1 , 




a a u=0.60 






* * U=0.64 _ 



















FIG. 2. Average relaxation time t av _ re i ax for an initial 
MShp configuration as a function of N for several values 
of U with a = 0.9. The upper solid line indicates the slope of 
the linear TV dependence for comparison. 
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FIG. 1. The entropy 5 as a function of U (upper graph) and 
the corresponding energy temperature relation (lower graph) . 
The equilibrium values of S and T vs. U is plotted as a 
solid line, both above (HP) and below (CP) the critical en- 
ergy U c = 3/4. These same quantities, evaluated for the 
metastable homogeneous phase MShp, are represented by 
dashed lines. The inset in the upper graph zooms in the en- 
ergy range where the phase transition occurs. The vertical 
dotted line indicates the critical energy U c - 
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